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Abstract 

We investigate baryon and charge transport in relativistic heavy-ion colh- 
sions, compare with Au + Au RHIC data at ^/snn = 0.2 TeV, and make 
predictions for net-proton rapidity distributions in central Pb + Pb coUisions 
at CERN LHC energies of ^s]^ = 2.8, 3.9, and 5.5 TeV. We use the gluon 
saturation model and put special emphasis on the midrapidity valley \y\ < 2. 
Net-kaon distributions are calculated and compared to BRAHMS Au + Au 
data at RHIC energies of y/s^jy = 0.2 TeV, and predicted for Pb + Pb at 
5.5 TeV. 

Key words: Relativistic heavy-ion collisions, Net-proton and net-kaon 
rapidity distributions. Saturation model. Predictions at LHC energies 
PACS: 24.85.+P, 25.75.-q, 25.75.Dw, 12.38.Mh 

The test of gluon saturation in relativistic heavy-ion collisions is an im- 
portant aim of the forthcoming Pb + Pb experiments at the LHC. At these 
energies gluons dominate the dynamical evolution of the system, which is 
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driven by a single hard scale, the saturation scale Qs ^ ^qcd [li- Whereas 
most theoretical investigations concentrate on charged-hadron production 



from inclusive gluon interactions 



3|, the valence-quark scattering off the 



gluon condensate as an observable in net-baryon distributions is expected 
to provide interesting new information on gluon saturation, and on geometric 
scaling 

Here the most promising effects arise at very forward angles, or corre- 
spondingly large values of the rapidity y ~ 5 — 8 at LHC energies of ^/snn = 
5.5 TeV for Pb + Pb, with a beam rapidity of 8.68. For symmetric systems, 
two symmetric fragmentation peaks are expected to be present in the net- 
baryon distributions at forward/backward rapidities. In particular, we have 
shown in that it is in principle possible to determine the growth of the 
saturation- scale exponent, A = d\nQs/dyh, with the beam rapidity yb from 
the position of the fragmentation peak in rapidity space. 

In the region of relatively large values of Feynman-x ^0.1, the valence- 
quark parton distribution in the projectile is well-known close to and below 
its maximum, and can hence be used to access the gluon distribution at 
small X in the other nucleus where saturation is expected to occur due to the 
competition of gluon recombination with the exponentially increasing gluon 



splitting 



mm. 

Whereas it is interesting to investigate these effects theoretically, the 
forthcoming LHC experiments with heavy-ion capability ALICE, CMS, TOTEM 
and ATLAS initially will not be able to detect identified baryons and an- 
tibaryons from central heavy-ion collisions at large values of rapidity in the 
fragmentation-peak region. In particular, the dedicated LHC heavy-ion ex- 
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periment ALICE will provide full particle identification for protons and an- 
tiprotons as well as kaons only in the central part of the net-baryon midra- 
pidity valley [9|. 

As a consequence, one first has to concentrate on the midrapidity region in 
order to compare with data. In this Letter we present our predictions for net- 
baryon rapidity distributions in central Pb + Pb collisions at LHC energies. 
The LHC physics program starts with center-of-mass proton-proton energies 
of 7 and possibly 10 TeV, the corresponding energies for Pb + Pb (scaling 
with Z/A) are ^/sNN = 2.76 and 3.94 TeV. We also present predictions for the 
highest attainable Pb + Pb energy of 5.52 TeV. Since experimental results 
will be available for net protons, we calculate these at the highest LHC energy 
in the midrapidity valley \y\ < 2 instead of net baryons, and also include a 
prediction for net kaons {K^ — K^) since these carry part of the valence 
quarks. 

The differential cross section for valence quark production in a high- 
energy nucleus-nucleus collision is calculated from 

xiq^{xi,Qf) ip{x2,PT) , (1) 



cPprdy (27r)2j9^ 

where pt is the transverse momentum of the produced quark, and y its 
rapidity. The longitudinal momentum fractions carried, respectively, by the 
valence quark in the projectile and the soft gluon in the target are Xi = 
Pt/ Vsexp{y) and X2 = pr / y/s exp{—y) . The factorization scale is usually 
set equal to the transverse momentum, Qj = pr ilO|. We have discussed the 
gluon distribution (p{x,pt) and details of the overall model in j4|. 

The contribution of valence quarks in the other beam nucleus is added 
incoherently by changing y — )■ —y. The valence quark distribution of a 



nucleus, = q — q, is given by the sum of valence quark distributions q^^^ 
of individual nucleons, q^ = Aq^^N, where A is the atomic mass number. 

Assuming that the rapidity distribution for net baryons is proportional to 
the valence-quark rapidity distribution up to a constant factor C, we obtain 
by integrating over pt, 

dN C f (fpT . ^ ^ . ^ 

^ = ^y ^x,g.(xi,Q,)v.(x„p^). (2) 

It turns out that this is indeed a good approximation at sufficiently high 
energy, in particular, when comparing to Au + Au data from RHIC, and we 
expect it to be valid at LHC as well |4|. 

The unintegrated gluon distribution is peaked at qr = Qs, or xi = 
exp (— r/2 + A), with the saturation momentum squared = A^^^Qlx2^, 
the saturation-scale exponent A, and the scalingvariable r = ln(s/(5o) ^ 
In ^1/3 _ 2(1 + \)y that we have introduced in ^. Here A is the nucleon 
number and Qq sets the dimension. The peak at g-r = Qs reflects the fact that 
most of the gluons sit at this value. Therefore, we expect dN/dy ~ Xiq{{xi)), 
with (xi) = {Q s) / ^ ex'^iy) . With X2 = xiexp{—2y) we can solve this equa- 
tion for yielding 



{xi) = { ^ ) exp 



'2Ta^ 



(3) 



In the region of small Xi and X2 corresponding to the midrapidity valley 
{y ~ 0) away from the peaks, the valence quark distribution behaves as 
, where the intercept A has been calculated in the saturation picture 



11] leading to 



with Cf = {Nl - l)/2Nc, Nc = 3. 
the old prehminary BRAHMS data in 
strong-couphng constant as — 0.3. 

Therefore, in the midrapidity valley Eq. ([2]) becomes 



he value of A had been fitted to 



ll|, with A ~ 0.47, leading to a 



1 dN MVSQo 
oc 



A dy 



/I 



A/(l+^) 



cosh 



1 + A 
2A -y 

2 + 



(5) 



which reduces to Eq. (80) in for the special case A = 0. The midrapidity 
values of the net-baryon or net-proton rapidity distributions at two different 
center-of-mass energies in the nucleon-nucleon system are related through 



dN 
dy 



s) 



A/i2+X) dN ^ 

-J) 



(6) 



We now use the analytical form for dN/dy = acosh.{by) (cf. Eq. ([5])) in a 
direct comparison with BRAHMS net-proton data in central Au + Au colli- 
sions at RHIC energies of ^/snn = 0.2 TeV |l2| through a x^— minimization 
of the two parameters a and b, where b = 2A(1 + A)/(2 + A). 

Our comparison with the BRAHMS Au + Au data in the midrapidity 
region is shown in Fig. [T]for net protons. The fit parameters are a = 6.79 ± 
0.59, b = 0.575 ± 0.116 as results of the x^— minimization per degree of 
freedom (8 data points - 2 free parameters; x^/dof = 0.028). With the energy 
dependence as expressed in Eq.Q, the grey band in the lower part of the 
figure shows our predictions for central Pb + Pb collisions at a LHC energy 
of ^/sNN = 2.76 TeV with A = 0.3 (upper bound), A = 0.2 (solid curve), and 
A = (lower bound). The mass-number dependence is very weak, and we 
neglect it in the discussion {Apb/AAu — 1.056). 



Note that for A = 0, we have A = b. Our value for b is slight 
than, but within our error bars compatible with the one fitted by 



arger 
. We 




Figure 1: The rapidity distribution of net protons in central (0%-5%) Au + Au collisions 
at RHIC energies of ^/snn = 0.2 TeV as measured by BRAHMS (black dots) is 
fitted with our theoretical formula using a x^— minimization to fix the parameters for the 
predictions at LHC energies. The data point at y = 2.9 is neglected in the minimization. 
The grey band in the lower part of the figure shows our predictions for central Pb + Pb 
collisions at LHC energy of y^s^^ — 2.76 TeV (corresponding to 7 TeV in p + p) with 
A — 0.3 (upper bound), A = 0.2 (solid curve), and A = (lower bound), using Eq. 
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Figure 2: The net-kaon rapidity distribution in central (0%-5%) Au + Au collisions at 
RHIC energies of y/s^N = 0.2 TeV as measured by BRAHMS |l2| (black dots) is fitted 
with our theoretical formula using a minimization. The grey band in the lower part 
of the figure shows our predictions for central Pb + Pb collisions at LHC energy of ^/snn 
= 2.76 TeV (corresponding to 7 TeV in p + p) with A = 0.3 (upper bound), A — 0.2 (solid 
curve), and A = (lower bound). 
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extract a value of A 0.575 for A = 0, and A ^ 0.509 for A = 0.3, leading 
to — 0.2 in both cases due to compensating effects. 

Our result for the midrapidity distributions should be compared directly 
to the forthcoming ALICE net-proton data in central Pb + Pb collisions. 
The predicted midrapidity value at ^/sNN = 2.76 TeV is dN/dy ~ 1.93, 
it depends only slightly on the saturation-scale exponent A and hence, one 
cannot expect to determine the value of A from midrapidity net-baryon data. 

For a determination of A from heavy-ion data at LHC energies, one there- 
fore has to rely on the forward rapidity region for net baryons [4| which is 
difficult to access experimentally, or on midrapidity distributions for pro- 
duced charged particles [3|. From the overall accuracy of our prediction 
regarding the absolute value at midrapidity, and the shape of the net-proton 
rapidity distribution, we will, however, be able to draw conclusions regarding 
the validity of the gluon saturation picture. 

In the comparison of our model calculations with RHIC net-baryon data 
in ^, medium and final-state effects turned out not to be important for 
rapidity distributions, although they may be visible in transverse momentum 
distributions. Due to baryon-number conservation, they could only lead to a 
redistribution in rapidity space, which is, however, not observed. This result 
differs qualitatively from the importance of medium effects in jet suppression 
at these energies. The expectation is that net-baryon distributions at LHC 
are also not visibly affected by medium and final-state effects. 

Like net baryons and net protons, net charged mesons such as kaons 
and pions carry part of the valence quarks, and can thus be treated on 
the same footing. In particular, these can be used as a cross-check for the 
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validity of our hypothesis that net-baryon and net charged-meson rapidity 
distributions essentially reflect the valence quark distributions, such that 
hadronization does not play a significant role. Here we study the net-kaon 
rapidity distribution since we don't have access to the full net charged-meson 
distribution. 

In Fig. [2] we show the result for the net-kaon rapidity distribution in cen- 
tral Au + Au collisions at RHIC energies of ^/snn = 0.2 TeV in comparison 
with BRAHMS data [l^ through a x^— minimization of the two parameters, 
a = 2.087 ± 0.173 and b = 0.535 ± 0.031. Here the result of the minimization 
is x^/dof = 0.540. Interestingly enough, the value of b for Ak is compatible 
with the one extracted for net protons. This indicates that the rapidity dis- 
tribution is primarily sensitive to the initial conditions of the collision, not 
to the hadronization process, since the slope is not depending on the species 
of the produced particles (protons or kaons). 

We had discussed the effect of fragmentation in our earlier work [4] by 
comparing calculations for net baryons with and without fragmentation func- 
tion, cf. Fig. 5 in that work. In the midrapidity region \y\ < 2, the effect at 
y/sNN = 0.2 TeV is clearly smaller than the size of the experimental error 
bars. We have therefore not discussed the effect explicitly in the present 
work. Again our extrapolation to central Pb + Pb at a LHC energy of 2.76 
TeV is shown by the band in the lower part of Fig. 2, with the solid curve 
for A = 0.2. Our predicted midrapidity value is dN/dy ^ 0.7. 

In Fig. [3] we display the energy dependence of our net-proton central Pb 
-|- Pb results near midrapidity for yj/vjv = 2.76, 3.94, and 5.52 TeV. At ?/ = 
the corresponding values of dN/dy are 1.9, 1.7, and 1.4. 
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Figure 3: Rapidity distributions of net protons in 0%-5% central Pb + Pb collisions at 
LHC energies of y^SNN = 2.76, 3.94, and 5.52 TeV. The theoretical distributions are shown 
for A = 0.2. 



A description for the net-proton rapidity distribution within the rela- 
tivistic diffusion model (which is not based on QCD, but on nonequihbrium- 



statistical physics) had been developed in 13|. There the predicted midra- 



pidity value for central Pb + Pb at the LHC energy of a/JJvtv = 5.52 TeV is 
dN/dy ~ 1 - 2.5 depending on the model parameters and hence, comparable 
to the present Q CD-based result. 

To emphasize how our midrapidity results are embedded into the over- 
all shape of the rapidity distribution for net protons (baryons) in central 
relativistic Pb -|- Pb collisions at LHC energies, we show the total rapidity 
density distribution functions for the BRAHMS Au + Au data at 0.2 TeV 
and for the three LHC energies ^/s^ = 2.76, 3.94, 5.52 TeV in Fig. H 

Here we have used for the mid-rapidity valley Eq. ([5]) (as in Figs. 1, 2) 
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Figure 4: Calculated rapidity distributions of net protons in 0%~5% central Pb + Pb 
collisions at LHC energies of y^SNN — 2.76, 3.94, 5.52 TeV. Our result for central Au + 
Au collisions at RHIC energies of 0.2 TeV is compared with BRAHMS data 12] in a x^ — 
minimization as in Fig. [TJ 
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matched at the point X2 = 0.01 with the parametrization (cf. Eqs. (7) and 



(8) in 15| ) of the valence quark distribution function which describes the 
larger rapidities. As is evident from the figure, the transition between the 
two regimes is fairly smooth. Both up- and down-quark parton distribution 
functions are considered with the appropriate weights. 

To conclude, we have presented predictions for net-proton and net-kaon 
rapidity distributions in central Pb + Pb collisions at LHC energies, with 
an emphasis on the midrapidity region where data will be available in the 
near future. We have set up and used a transparent QCD-based model and 
well-established parton distribution functions in the context of saturation 
physics, and we expect that our predictions at midrapidity will turn out to 
be reliable. 

We have shown that hadronization does not influence the slope of the 
net-hadron rapidity distributions since net-proton and net-kaon rapidity dis- 
tributions are related through a constant factor. Hence, net-baryon and 
net-charge transport provide a powerful tool to investigate initial-state dy- 
namics in heavy-ion collisions. Finally, we have extracted a value for the 
strong-coupling constant ~ 0.2 both from net-proton and net-kaon rapid- 
ity distributions in Au + Au at RHIC energies. 

This work has been supported by the ExtreMe Matter Institute EMMI. 
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